ON POSITIVE CLASSICAL AND FREE STABLE LAWS 



N. Demnfl 



Abstract. We derive the representative Bernstein measure of the density of 
, < a < 1, where X a is a positive stable random variable, as 
a Fox-H function. When 1 — a = l/j for some integer j > 2, the Fox H- 
function reduces to a Meijer G-function so that the Kanter's random variable 
(see below) is closely related to a product of (j — 1) independent Beta random 
variables. When a tends to 0, the Bernstein measure becomes degenerate 
thereby agrees with Cressie's result for the asymptotic behaviour of stable 
distributions for small values of a. Coming to free probability, our result 
makes more explicit that of Biane on the density of its free analog. The paper 
is closed with analytic arguments explaining the occurence of the Kanter's 
random variable in both the classical and the free settings. 



1. Motivation: Kanter's random variable 

The study of stable random variables has known a considerable growth during 
approximately the last five decades, since they are shown to be an efficient model 
for various phenomena occuring in biology, quantum physics, market finance f [16]). 
Unfortunately, their densities has no closed formulas except in few cases. Consider 
for instance, the stable random variable X a of index < a < 1 and of asymmetry 
parameter equal one, known also as the positive stable variable since it is supported 
by (0, oo) ([[5J, p. 50). Its Laplace transform is given by e - *" ,t > and its density 
admits the following expansion ((2.41) p.54 in [9], (14.31) p.88 in [12], see also [17] ) : 

M) fe x , v 1 



X ( — 1) 

fe=l 

The positivity of this series is not trivial and follows for instance from the integral 
representation displayed p. 74 in [16 . Amazing and not trivial as well is the com- 
plete monotonicity of the density of (X a )~ r , r > a/(l — a) (in particular its infinite 
divisibility, [12]) which follows from the following representation ([H]): there ex- 
ists a standard exponential random variable L and an infinitely divisible random 
variable V a ^ r , both variables being independent and 



e v a 



Recall that the proof of this fact relies on a suitable integral representation of the 
Gamma function and that the random variable V a ,r is characterized by its Laplace 
transform. For r = a /{I — a), the above representation is more precise and traces 
back implicitely to Chernin and Ibragimov's paper [4] (see also Ch.II in [9]). In 
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fact ([10, p.703) 

(Y W(l-«) — 

where U is a uniform random variable on (0, 7r) and 

a a {u) = [b a {v)]W- a \b(u) = [c a (u)] a [c^ a {u)]^- a \c a {u) := 

sm(u) 

We shall refer to a a (U) as the Kanter's random variable. From the injectivity of 
the Laplace transform, one gets 

V a ,i- a = log[a a (U)]. 

Mysteriously, the function a a is intimately related to the density of the free stable 
law of index < a < 1 and asymmetry parameter p = 1 derived p. 1053 in [T], that 
we shall call shortly the positive free stable law. More precisely, if U is defined on 
some probability space (f2,jF,P), then 

P{a a (U) <u) = -tC 1 ^), u > a a {0) = (1 - a)a Q/(1 ~ a) 

since a a is strictly increasing ([10] p. 704), while the density of the positive free 
stable distribution takes the form ([I]): 

1 sm(a^ a {x))sm(aa^ a {x)) 



(1) 



sin((l - a)a 1 i Q (x)) 



for x > ai_ Q (0). That is why we thought it is reasonable to start with seeking 
an expression for the density of a a {U) as a special function and more generally for 
that of e Va r . The easiest way to proceed is to invert the Mellin transform of the 
density of e~ Va - r given by ([H p.292): 

T{rs/a + 1) 



E(e S r( a + i)r(«r + i)' 

However, the Mellin's inversion formula 



U(s) > -a/r > -(1 - a). 



(2) 



1 



2iny 



-\-ioo 



T(rs/a + 1) 

r(s + i)r(sr + i 



■y s ds 7 



applies for any r > a/(l — a) (we take 5ft(s) = as a path of integration) and does 
not when r = a/(l — a). The latter fact is seen from the following estimate of 
\T(a + ib)\, a, b e K, |6| — > oo displayed in formula 1.24. p. 4 from [8]: 



(3) 



\T(a + ib)\ - V%ir\b\ 



a-l/2 -o-7r|6|/2 



\b\ 



oo 



so that the absolute value of the integrand in © is equivalent to |s| -1 / 2 (up to a 
constant). Note however that when r > a/(l — a), the integral displayed in ([2} is 
nothing else but a so-called Fox H- function ([5] p. 3, see below): 



1 



+ioo 



T(rs/a+ 1) 



1 



-y-* d s = -H^ 



2inyJ_ loo r(s + l)r(sr + l)* y 
so that the density of e Va - T reads 



,.|(bi),(i,r) 

y\(l,r/a) 



1 



-H. 



1.0 

2,1 



1 



1(1, r/a) 



y > 0, r > 
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The main result of this paper states then that in the pathological case r — a/(l — a), 
the density of a a (U) is a H- function too represented through a path integral re- 
sembling to @ but with a different contour (Hankel path). When a approaches 
zero, this density tends to zero pointwisely, a fact that agrees with Cressie's re- 
sult implying that the Kanter's random variable becomes degenerate in the limit 
a — > + . When 1 — a = 1/j for some integer j > 2, the derived H- function reduces 
to the so-called Meijer G-function so that a a (U) is closely related to a product of 
suitably chosen j — 1 independent Beta random variables. The latter fact resem- 
bles to Williams respresentation of positive stable variables ([H])- Coming to free 
probability theory, the distribution function of a a (U) is computed yielding an ex- 
plicit expression for the density of a positive free stable law. Using the free analog 
of Zolotarev's duality ([I]), a similar result holds for a free stable law with index 
1 < a < 2 and asymmetry parameter p — 0. The paper is closed with supplying 
analytic arguments that explain the occurence of a a in both probabilistic settings. 
It turns out that it is resumed in the inverse formula for the Fourier transform of 
the density of X a together with deformations of paths of integration. Nevertheless, 
we think that a group theoretical argument exists but seems to be very hidden. 
For sake of completeness, some basic needed facts on the H-function are collected 
in the next section. A good reference to them is the monograph [3]. 



2. On the Fox H-function 
The Fox H-function is defined as a Mellin-Barnes integral: 

i(oi,Ai)i<i<j 



TT7n,n 



0{s)z- s ds 



L 



where 1 < m < q, < n < p, a t , b, e K, A i ,B l > 0, 



©00 - rW 



Ut m+ i r(i - h - Be) nf=„ + i r(a, + A iS ) ' 

the principal determination of the power function is taken (though it is not necessary 
here) , an empty product being equal one and L is a suitable contour separating the 
poles of both products of the meromorphic (Gamma) functions in the numerator 
of 9 ([8] p. 3, 4). The choice of L and the domain of convergence of the Mellin- 
Barnes integral defining the H-function for each contour depend on the parameters 
ai,bi, Ai, Bi. Eight cases are discussed in [S] p. 4. and we will only make use of five 
of them in the sequel that we shall refer to whenever needed: 

i) q> 1, 

q p 
i=l i=l 

then the H-function exists in the punctured complex plane. 

ii) p > = 0, the H-function exists in the domain \z\ > /3 where 

■ ■ II i 1 IK 

i=l i=l 

and L = L +OQ described in [5] p. 3. 
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iii) Let fi be the real number defined by 

romp q 
i—1 i—1 i—n+1 i—m+1 

Then the Fox H-fvmction exists for all z =^ such that | arg(z)| < 7rfi/2, f2 > 
and L is the infinite semi-circle 7 — ioo, 7 + ioo for a suitable real 7. 

iv) When q > 1, /x = 0, the H-function exists in the open disc |z| < /3 and 
L = L-00 is described in [8] p. 3. 

v) Let 

q p 

p-q 



2 

If ji = 0, 5 < — 1, then the H-function is well-defined for complex numbers 
lying on the circle \z\ = /3. 

The most crucial fact when looking at these definitions is that whenever more 
than one definition make sense, they lead to the same value of the integral. In 
fact, the integral is evaluated by means of Cauchy's residue Theorem applied to 
truncated contours together with asymptotic estimates at infinity, and the contour 
of integration can be deformed into another one provided that the deformation 
is admissible ([6], [14], see also CH.V. in [7 for a similar discussion on Meijer's 
G- function which fits the Fox H-function when Ai = Bj = 1, 1 < i < p, 1 < j < q). 

Remark 1. In the sequel, the complex variable z will take only strictly positive 
values thereby s 1— > z~ s is defined through the real-valued logarithm function. 

3. The density of a a (U) 

According to Kanter's representation, there exists a density h a supported by 
a Q (0,7r) = (a a (0),oo) such that 

^ V ' ■ sin(/fc7ra)r(l + ka)x k ^- a) = / e~ xy yh a {y)dy . 



kl 

k=l 

The main result of this paper is stated as: 

Proposition 1. The density of the Kanter's variable is expressed as: 



h a {y) = ~ H lf 



1 



(l,a/(l-a)),(l,l) 



2(1 - a) gl|0 

~ H 2-l 



1 



,(l,2a),(l,2(l-a)) 



y 2(l-a) 1(1,2) 



for y G [a a (0),oo[. 

Proof. We shall start from the density of (X a )~ a which has the expansion 

1 00 (--[) k 

V sin(fc7ra)r(l + ka)x k . 

airx f-^ k\ 
fc=i 

Using the mirror formula satisfied by the Euler's Gamma function ([7 p. 3), 

r(l + ka)T(-ka) = - I 

the density of X~ a is transformed to 

1 (_!)*= x k 

ax kl T(-ka) 
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and the infinite summation may start from k = since z = is a pole of the Gamma 
function. According to A. 6 p. 218 in [8], this is expressed through the H-function 
as follows 0: 



■E 



- — sin(fc7ra)r(l + ka)x k 
k\ 



-H 



i,i 



.2 1 (0,2a) 
1(0,2) 



for all x > (note that the H-function in the RHS exists for all z ^ in the complex 
plane by i) since q > 1, fi = 2 — 2a > 0). Therefore, the density of (A Q ) _Q ^ 1_a ) 
reads 

" 2(l-a)|(0,2a)" 



2(1 -a) „i, 
"-"1,1 



cv;r 



According to hi), one has 



H 



i.i 



J2(l-a)|(0,2o) 
1(0,2) 



1 

2i7r 



1(0,2) 



7+100 r(2a) da 
7 „ loo T(2as) a;2(i-«)* 



where 3?(s) = 7 > 0. With the help of the Gamma integral 

2(1 -a) n +lo ° r(2s) 

/ T -ioo r(2a S )r(2 + 2(l-a) S ) ./„ 



2(1 - a) 
ax 2 M 



„2(l-a)|(l,2a) 
1(1,2) 



2m7r 
2(1 -a) ri+' x - 



2l7T 



T(2s + 1) 



7 — ZOO 



r(2as + l)r(2 + 2(1 - a)s) 



e- xy y 2{1 - a)s+1 dyds 



e- x yy 2 ^-^ s+1 dyds. 



In order to change the order of integration, one needs to check that 
+oc |r(2is + 2 7 + l)| 



xy, .2(1— a)7 



T{2ias + 2a7 + l)r(2(l - a) a + 2(1 - a) 7 + 2; 
Performing an integration with respect to y, it amounts to check 

+°° |r^ s + 2 7 + i)| 

.co \T(2ias + 2aj + l)r(2i(l - a)s + 2(1 - a)j + 2)| 



dyds < 00. 



ds < 



which follows again from the estimate of \T(a + ib)\, \b\ — > 00 (see ©). Thus, 
Fubini's Theorem yields 



J yh a {y)dy = 2(1 - a) / xe x 



(4) 



7+ioo 

7 — ZOO 

1.0 



T(2S + 1) 



r(2as + l)r(2 + 2(l-a)s) 



2(l-a)s+l 



dsdy 



/>oo 

2(1 -a) xe-^yH 1 ^ 
Jo 



1 



2(l-a) 1(1,2) 



i(l,2a),(2,2(l-a)) 



dy. 



Note that the H-function displayed in ((3]) may be represented for y > (1— a)a a l 
as a path integral through the loop L_oo according to iv), it vanishes for < y < 
(1 — a)a Q ^ 1_Q ' since according to ii) the contour is a loop £+00 that contains no 
pole of the meromorphic function s T(l + 2s), and it takes a finite value at 
z = (1 — a)a Q /^ 1_Q ^ regarding v) since /1 = 0,6 = —3/2 < —1 (in particular, it 
vanishes there by continuity). Hence, an integration by parts shows that the density 



! It is a special kind of what V. M. Zolotarev called incomplete hypergeometric function in I17| . 
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of the Kanter's random variable is 



h a (y) 



2(1 -a) d 
V dy 
2(1 -a) d 



i 



I (1,2a), (2, 2(l-a!)) 



2(1-0) 1(1,2) 

r(2« + i) 



V dy Jr T(2as + l)r(2 + 2(1 - a) a) 



2(l-«) S +l ds 



2(1 - a) 



r(2s + 1) 



i _ oo T(2as + l)r(2(l - a)s + 1) 



2(l-a) Sds 



2(1 i,o 

—y— H ^ 
1 



i 



,(1,200,(1,2(1-00) 



2(l-o) 1(1,2) 



y 

r( a /(l-a) + l) _irrl,0 

r( as /(i-a) + i)r( s + i) y 2 ^ 



,(l,a/(l-a)),(l,l) 
l(l,l/(l-o)) 



where L'^^ is the image of £_oo under the map s n- 2(1 — a)s. The density /i a is 
well defined for y > (1 — a)a a ^ 1_Q ' ) by iv) and vanishes for < y < (1 — o^a ^ 1-0 ) 
according to ii) (see also 1.33 p. 6 in [5]). ■ 

4. Special indices 

4.1. The 1/2-stable case. When a = 1/2, the density of the positive stable ran- 
dom variable reads ([16] p. 66) 

1 T -3/2 -l/(4x)i 

According to Kanter's representation, the image of the above density under the 
map x i — ^ 1/x, that is 



l - x /4i 1 r^t 



x/4 



should be the Laplace transform of y t— > yhi/ 2 (y) which takes the form (after the 
use of the mirror formula) : 



yh 1/2 (y) = ^2 E sin2 t( fc + W^lfr + !)/ 2 - 



(-1)* 



fc=0 



k\ y( k + 1 )/ 2 



, y > 1/4. 



Note that x H> a; : / 2 is completely monotone since x H> x 1 / 2 is a Bernstein function, 
and that the representing measure is easily computed using the Gamma integral as 

dy 

\y>0}- 



Now, since sin(n7r) = for any integer n, then 

1 

2^2 



, , , 1 ^sin 2 [(2fc + l)V2]r 2 [(2fc + l)/2] 1 



A:=0 



1 °° 
2, ,1/2 



(2*0 

r 2 (fc + i/2) i 

2vrV/2 ^ r(2fe + l) 
for y > 1/4. Using Legendre's duplication formula ([7] p. 5) 

v / 7rr(2fc + 1) = 2 2fc r(/c + i/2)r(fc + 1), 



y(2fc+l)/2 
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one gets 



yhi/ 2 (y) 



1 X 
-,,1/2 



(1/2). 



^ fco fc! (4y)fc 27ryi/2 ^ _ 1/(4y)) 
1 1 

y > 1/4, 



which is nothing but 



2^7^ 1{y>0} * ,5l/4(d2/) 



where * is the classical convolution of positive measures. 

4.2. Relation to product of independent Beta variables. Let 

1 3 



1 — a = — , a 
3 



3-1 



for some integer j > 2, then 
E 1 



K-i/^)] 5 



r(j S + 1) 



r( a + i)r((j - i) s + 1) 



»(*) > -i/i- 



The multiplication Theorem (generalization of Legendre's duplication formula, [7 
p.4) 



T(js + 1) 



1/2 rj r 



fc=i 



s + - 
3 



yields 



E 



3-1 



Z TT 

2tF V J - 1 (3 - 1) (J - 1)S H r(a + k/(j - 1)) ' 



r (s + k/j) 



Now, let /3i, /?2, • • • , /3j-i be independent Beta random variables of the first kind 
such that the density of /3k reads: 

where B(-, •) stands for the Beta function. According to formula (4.11) p. 122 in [8], 

r(a + fc/j) 
11 T(k/j + 1) 11 r(a + k/(j - 1)) 

1 



/i-i 



r(j - 1) 



2tt 
1 



n 



r(a + fc/j) 



j (i-i)O'- 1 ) ry) a = i r( fl + fc/(j - 1)) 
'' r(s + k/j) 



i 



- 1)0-1) 11 



again by the multiplication Theorem. It follows that 



\ T(s + k/(j - 1)) 
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When s ^ 0, the last equality may be written as 

w s_x P(a 1 _ 1/i (U') < l/u)du = jE 



(j- 1)0-D 



s-1 



where V is uniformly distributed in (0, 1) and is independent from the /?fcS. 

Remark 2. The relation to product of independent Beta variables is not surprising 
since the H-function involved in the density ofhi_ 1 /j reduces, when 1 — a= 1/j to 
a Meijer's G-function ({5] p. 123). It somewhat matches with Williams result f\15\ ) 
if one keeps in mind the Gamma-Beta algebra. In fact, 1/X a is distributed, when 
a = l/j,j>2asa product of j — 1 Gamma random variables 7fc, 1 < k < j — 1 
with densities 

~'"L{t>0}> l<fc<j-l- 



nk/j) 



t*/i-i e -*l. 



Moreover, since 



K(0)} 



{l-a)/a 



f in a Williams result implies 



fix 



ai- a (6), 



■Ik 



Xi/j a i-i/j(U) 

Note also that, for a = l/(jd),j > 2d,d > 1, a more general representation of 
positive stable laws is due to T. Simon (\Vd\) and involves independent Gamma and 
Beta variables. 

5. Free positive stable distribution 

Free stable distributions may be defined in a similar way as in the classical 
setting, when substituting the classical convolution of probability measures by 
Voiculescu's free convolution (see [1]). As mentioned in the introduction, the in- 
verse function (in the composition's sense) of a\- a provides an explicit expression 
of the density of a stable law of index < a < 1 and asymmetry parameter p = 1 
in the free probability setting. More precisely, if V is a uniform random variable 
on (0, 1) defined on some probability space (f), J?,P), then 



for x > a(l 



-a^ a (x)=V{ ai - a (TrV) < x) 

IT 

a )(i-<*)/«. Wnen a = i/2, it simplifies to 

dy = 1 r~ 1/A 



"1/2^ = 



= 2 



/ 4 y\/y- 1/4 
2dy 



dy 



(y + l/4)y/y 



1 + Ay 2 

Using the trigonometric identity 

sin(0) = 2 



= 2 arctan 



tan(6»/2) 



1 + tan 2 (6»/2) 
one recovers the density given by Biane ([T]) 



sin[a 1/ 1 2 (x)] 



tan(a 1/ 1 2 (a;)/2) 



1 



ttx 1 + taxi 2 {a~Ux) 1 2) ttx! + i(x - 1/4) 



2ttx 2 
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For general a £ (0, 1), we use the integral representation 



H. 



1.0 



1 



,(l,2a),(l,2(l-a)) 



2a 1(1,2) 



T(2S+ 1) 



y 2as ds. 



2inJ L _ Y(2as + l)r(l + 2(1 - a)s) 



Now, assume for a while that Fubini's Theorem applies so that 



r(2s + i) 



I - 1 ( \ = — 

^1-aW 2ijT T ^ 2as + 1)r(1 + 2 ^ _ ^ 

r(2,s) 



y 2as - x dyds 



2i7r Jr T(2as + 1)T(1 + 2(1 - a) a) 



(l_ Q )(l-»)/c 

(a 2a {l- a)^ L - a >) 6 



ds. 



Since the poles of s i— > T(2s) are s — —k/2, k £ N and simple, and since the pole 
s = lies outside -L_oo, then Cauchy's Residue Theorem yields 



1 

2i7T 



T(2s) 



x 2as ds = Hff 



T(2as + l)r(l + 2(1 - a)s) 



, (1,2a) ,(1,2(1-0:)) 



x 2a 1(0,2) 



-Res(x, 0) 



where Res(:r, 0) is the residue of the meromorphic function 



r(2*) 



„2as 



r(2as + l)r(l + 2(1 - a)s) 
at s = for fixed parameter x > a{\ — a)( 1 ~ a ^ a . This residue is easily computed 



lim 



sr(2s) 



>a 2a (l-a) 2 ^- a K 



>or(2as + l)r(l + 2(l-a)s) 
Since the last H-function is defined on the whole real line, in particular at a(l 
q,^ 1 -")/" by th e virtue of v) since p = and 5 < —3/2 < —1, then 



-a^ a (x)^2H^ 



1 |(l,2a),(l,2(l-q)) 
2,2a 1(0,2) 



1 



i 2a( 1 _ a )2(l-a) 1(0,2) 



,(l,2a),(l,2(l-a)) 



Regarding the continuity of the H-function at ail — a)' 1 Q )/ Q 5 it vanishes there 
(since it vanishes for x < a(l — a)( 1 ~ a ^ a ) and one gets 



~a^ a (x) = 2H# 

IT 



1 



,(l,2a),(l,2(l-a)) 



r 2a 1(0,2) 



,x>a(l-a) (1 ~ a)/Q 



and zero otherwise. Finally, in order to fit in the density of the free positive stable 
law displayed in ([1}, we can rewrite the last H-function as 



a* 2 '' 



l,(l,l),(l,(l-a)/a) 
3-1(0,1/0) 



1 



1 — a 



o-i,o 

-"2,1 



1 



,(l,a/(l-a),(l,l) 



c a/(l-a) l(0,l/(l-a)) 



Note that the duality relation given in [T] allows to derive the density of a free 
stable distribution of index 1 < a < 2 and asymmetry coefficient p — (the case of 
an asymmetry parameter p = 1 is obtained by the simple variable change x <-¥ —x). 
Coming back to the validity of Fubini's Theorem, one has to prove, after integrating 
with respect to the variable y, the convergence of 



|r(2s + 1) 



L _ \T{2as + l)r(l + 2(1 - a)s) 



c 2aSR W _ (a 2a (1 _ a) 2(l-a) ) K W 



\ds\ 



This is easily cheked from the estimate of T(a + ib) for large a (formula 1.23 p. 4. 
in!): 

|r(o + *6)| ~ V2^\a\ a - 1/2 e- a - a( ' 1 - asa W/ 2 , \a\ -> oo 
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which shows that the integrand is equivalent to 



C e -(alog(a) + (l-a)log(l-a))Sft( S )^ l^/xl _^ ^ 



|5R(s)| 3 / 2 

for some constant C = C(S*(s)) depending only on 3(s). 

5.1. Behaviour for small indices: Cressie's result and its free analog. In 

this section, we investigate the limiting behaviour of h a as a — ¥ 0. In the clas- 
sical setting, the density tends to zero pointwisely therefore agrees somehow with 
Cressie's result we recall below. The latter implies that a a (U) becomes degener- 
ate when a approaches zero. In the free setting, the limiting density is an infinite 
positive measure whose image under the map x M- 1/x is the Haar measure on the 
compact interval [0, 1]. Start with the series expansion of h a is given by (A. 6 p. 
218 in [§]) 

1 (-l) fc 1 

MW -(!-«) Z-^r(l-a-ka)r(a-k{l-a)) k\ y (fc+i)(i-a)+i 
for y > (1 — a)a a /^ a \ The mirror formula 

i> - Mi - «))r((i - a )(k + 1)) = — - — ^ 

sm((l — a)(k + l)n) 

transforms h a to 

1 - a ^ sin[(l - a){k + l)7r]r[(l - a){k + 1)] (-l) fc 1 
n ho T(l-a-ka) k\ y (k+x)(i-a)+i ■ 

The defining term of the last series converges as a — > to 

(-l) fc sin((fc + l)7r) 



yk+2 



= 



for any integer fc > while the support ((1 — aja"/' 1 a \oo) of h a decreases to 
[1, oo) as a — > 0. Since 

r[(l-a)(fc + l)] 
r(l -a- fca)fc! 

as a — > and it defines a continuous function on the variable a, then Lebesgue's 
convergence Theorem ensures the convergence of h a (jj), for fixed y > 1, to zero as 
the positive index a does. According to [5], (see also [3]), a stable variable with 
shape parameter c = 1 and zero shift 7 = satisfies 

so that a a (U) converges in distribution to the Dirac measure 5\. For a positive free 
random variable, we start from the image of its density function under the map 
x n- x a : 

a^l a (x^) sin( a ar-j£ 1/Q )) sin(a^ a (x^)) 



a a 
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for x > [ai_ Q (0)] a = a a {\ - a) 1 "" > 1/2. Then, we expand 



Z - ltl 2.1 



^ T(l 



1 |(l,2a),(l,2(l-cQ) 
x 2a 1(0,2) 

1 



k=0 



(l-ka)T(l-k(l-a)) k\ x k ' 



With the help of the mirror formula, 



r(l - fc(l - a))r((l - a)fc) = — — — , k>l, 

sm(/e(l — a)ir) 



one easily derives 

V- q (^ 1/q )-i+ 

a— >0 7T 



^ oo 

— lim y 

IT a-s-0 ^— ' 



sin(fc(l - a)7r)r(fe(l - a)) (-l) fc 1 
r(l — fca) kl x k 

for fixed a; > 1. Since a a (l — ->lasa^0, then 

a r- Q (^ 1/Q )sin(aa^ ct (^ 1/Q )) sm(a^ a (x^ a )) 1 



fc=i 

\ 1 — a 



lim 

a->0 



hia^ 1 ^) sin((l-a)ari Q (^)) 



for any x > 1. Note that Scheffe's Lemma does not apply as in the classical setting, 
nevertheless the image of 

dx 

x {x>1} 

under the map x i— > is the Haar measure on [0, 1], namely: 

dec 

— 1{0<z<1}- 



6. On the function a a 

As stated in the introductory part, the occurence of a a in both probabilistic 
settings is resumed in the integral representation of the density of X a together 
with deformations of integration paths. Start with the inversion formula for the 
Fourier transform to see that the density of X a reads ([9] p. 49): 

1 f°° 1 



31 / e- itx e- tae dt = -« / e ltx e- {lt) " dt. 



o 



The path of integration may be deformed into the positive half imaginary-line 
yielding ([H]) 



1 f°° 1 f°° 

-« / e ttx e- (lt)a dt = 3 / e- tx e~ eia 

n Jo * Jo 



dt 



-$S I e- tx e~ e *"**"<ft. 



Performing the variable change t = 

2,1/(1-00 

the last integral transforms to 

-l/fl-a) /-oo _l/(l-a) /-oo 
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where <j)i- a {u) := e l7ra u a . Following [4] (see also Ch.II [9]), the path of integration 
may be taken as the curve 

C a :={u = re t9 ,6e [-n, 0], S>(u + </>i- a (u)) = 0} 



u = re* 9 , 6> G [-7T,0],r = 
u = -re l6 ,e G [Q,7r],r = 



sin(a(7r — 0)) 



l/(l-a)- 



sin(ai 



sin 

1 V(l-«)' 



sm f 



Then for u 



G C a , 



3?(u + ^i_ Q (u)) 



■r cos 6 + r a cos(a6) 
*[-r 1_a cra0 + cos(a0)] 

a/(l-a) 



sin(a6 l ) 
sin 



sin((l - a)0) 
sin(0) 



and a a (6) shows up. Now, consider F a := 1/G a where G a is the Stieltjes transform 
of the positive free stable distribution ([I]). Then, F a is a one-to-one correspondence 
from the upper half-plane C onto its image. Let F~ l be the reciprocal function of 
F a , then (ITJ.i) 

F- 1 (z) = Q (z) + z 

where z belongs to the domain Q(z + 4> a (z)) > and $s(z) is positive and sufficiently 
large. The function <f) a is known as the Voiculescu transform of the positive free 
stable distribution. But since Stieltjes inversion formula involves the imaginary 
part of G a when z G C + approaches the real line, one rather focuses on the domain 

n a := {z, ^{l/z + (j) a {l/z) > 0} 

for z G C~ near the origin of the complex plane. This is a Jordan domain, that is 
delimited by a Jordan curve, whose boundary dQ a is obviously 

l/a ' 

dn a = { n"\tl ■ [-7r,0],r= ' " "' 



sin(7r + 6) 



sin((l - a)(Tr + 6)) 



G [0,tt], r 



sin^ 



sin((l - a)6)) 



l/a ' 



This is the image of Ci_ Q under the inversion map z t— > l/z. By Carathcodory's 
extension Theorem, G a extends to a homeomorphism from the real line onto dQ a 
so that every x lying in the support of the free positive stable variable is given by 

x = G- 1 (z), zedn a 

= 3t{l/z + </> a (l/z)). 

The same computations performed in the classical setting lead to the function ai- a . 
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